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Using a single functional form which is able to represent several different classes of 
statistical distributions, we introduce a preliminary study of the ferromagnetic Ising model 
on the cubic lattices under the influence of non-Gaussian local external magnetic field. 
Specifically, depending on the value of the tail parameter, r (r < 3), we assign a quenched 
random field that can be platykurtic (sub-Gaussian) or leptokurtic (fat-tailed) form. For 
r < 5/3, such distributions have finite standard deviation and they are either the Student- 
t (1 < t < 5/3) or the r-distribution (r < 1) extended to all plausible real degrees of 
freedom with the Gaussian being retrieved in the limit r — > 1. Otherwise, the distribution 
has got the same asymptotic power-law behaviour as the a-stable Levy distribution with 
a = (3 — t)/(t — 1). The uniform distribution is achieved in the limit r — > oo. Our results 
purport the existence of ferromagnetic order at finite temperatures for all the studied values 
of r with some mean-field predictions surviving in the three-dimensional case. 
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I. Introduction and motivation 

Disordered magnetic systems have long ceased being a "mere" theoretical exercise 
to find a limelight place in the class of problems which enables the reasoning about 
the equilibrium and out-of-equilibrium behaviour of several compounds that can be 
experimentally studied. Accordingly, the Sherrington-Kirkpatrick spin-glass model [1] 
and the random field Ising model (RFIM) [2] are quintessential systems that have soon 
found a correspondence with condensed matter and materials science case-studies [3]. 
For uncertain reasons, perhaps related to the nature of the replica method and the 
took for granted universality of the Gaussian, most of the studies have been devoted 
to the analysis of systems with n-modal or n-Gaussian disorder in either the coupling 
constant or the local external magnetic field [31 O El El El El HD1 EH]- However, the 
last decade has been particularly prolific in examples of the ubiquity of non-Gaussian 
distributions in natural systems and man-made phenomena [121 d3] • Therefore, it was 
completely upheld the introduction of generalized disordered magnetic systems in which 
fat tailed and compact distributions in the local magnetic field [Hj and fat tailed (o> 
stable Levy) distributions in the coupling constant were considered [151 EH El EE] • As 
a matter of fact, earlier experimental studies on organic charge-transfer compounds 
like N-methyl-phenazium tetra-cyanoquinodimethanide (NMP-TCNQ), quinolinium- 
(TCNQ)2, acridinium-(TCNQ)2 and phenazine-TCNQ suggested the experimental 
validity of non-Gaussian (fat-tailed and uniform) distributions. [Tj5] Furthermore, in 
alternative quantities and systems the importance of non-Gaussian distributions has 
been also demonstrated. For instance, in Refs. [121 E] , it is shown that for the a-stable 
Levy spin-glass the local fields are not Gaussian distributed. 

Concerning the non-Gaussian RFIM, a recent work [13] analyzed the case 
of distributions with power-law asymptotic behaviour, namely the generalized r- 
distribution and the Student-i [20] in the limit of infinite-range interactions. In this case, 
interesting results were found, namely the appearance of a change of the concavity of the 
critical transition line for r > 2 and a corresponding divergence of the zero-temperature 
free energy per spin. Despite the fact that mean-field surveys are unquestionably 
important for a sketchy description, most of the systems are not experimentally studied 
nor naturally exist at a dimension equal or higher than the upper critical dimension. 
Consequently, you should turn our attention to systems with achievable dimensionality 
and worth of study, i.e., exhibiting critical behaviour. In this particular case, it is 
known that the RFIM cannot undergo a discontinuous phase transition for d < 2 
because of the existence of a single minimum energy state [21], whereas Landau-like 
energetic arguments purport the impossibility of a ordered phase at 2D [2], i.e., no 
phase transition at all. Alternatively, the field-theory approach based on the existence 
of the dimensional reduction does not hold for the RFIM [22] , thus leaving the existence 
of critical behaviour for 3D systems under dispute. 

In order to better understand this problem, we herein report preliminary results on 
the study of the critical frontier of the RFIM on cubic lattices with short-range (nearest- 
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neighbour) interactions within the same quenched disorder conditions of Ref. [14J. To 
that, we performed Monte Carlo (MC) simulations on lattices with sizes in the standard 
range. The figures obtained point at the existence of ferromagnetic order at finite 
temperatures for a wide range of the parameters of the random-field distribution and 
that in some cases continuous phase transitions exist with the critical temperature being 
defined by both the tail and the width of the distribution. 



II. Model and Monte Carlo Simulation 

We have considered the following Hamiltonian on a cubic lattice of linear size L: 

n = -jJ2 s ^-J2 H * s ^ W 

<i,j> i 

where Si = ±1 and the sum Yl<ij> being performed over the nearest-neighbour (NN) 
pairs of spins. J is the ferromagnetic interaction constant between NN and the random 
fields {Hi} are quenched variables ruled by a PDF that is defined by a parameter r 
(generic degree of freedom) . For r < 1 , 



1 - r 



r 



5-3r 



P^Hi) = W— B T "" ' [1 - B T (1 - r)Hi\ — , (2) 



7T 



(with \H\ < [B T (1 — t)J ) which is the generalized r-distribution, and for r > 1, we 
have 



-B T ^^ [l - B T {\ - t)H?}^, (3) 



3-T 



_2(r-l) 

which is the generalized Student-t distribution also widely named g-Gaussian within 
statistical mechanics [13]. By generalized we mean that the degrees of freedom, m and 
n, of t- and r-distributions are extended to the entire domain of feasible real values 
according to the relations r = (m + 3) / (m + 1) [m > 0] and r = (n — 4) / (n — 2) 
[n > 2], respectively. In Eqs. fl2J) and (EJ, T[.] is the Gamma function and B T is given by 

Br = \ 2 , (4) 

(3 - r) a* 

where a T is the width of the PDF. For r < 5/3 the width and the standard deviation, 
a, are naturally related by, 

(5 - 3r) a 2 = (3 - r) a 2 T . (5) 

In Fig. [TJ we show the PDFs 02]) and ([3]) for some values of r, in the linear-linear 
(left panel) and log- log scales (right panel). We can see the power-law behaviour for 
r > 1. 

For the sake of simplicity, we shall consider J = 1 hereinafter. We studied systems 
of L = 6, 8, 10, 12 and 16, with periodic boundary conditions and a random initial 
configuration of the spins by means of Monte Carlo simulations. The algorithm is 
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Figure 1: (Color online) Random- field probability distributions for some values of the 
parameter r (from bottom to top: r = 5/2, 3/2, 1, 1/2 and —20), in the normal (left 
panel) and log- log scale (right panel). We have used oy = 1 in all cases. 

as follows: a configuration of the random field {H} is generated at the beginning of 
the simulation according to the corresponding probability Pi )S (H), and remains frozen 
during the whole dynamics; then, every lattice site is visited and a spin flip occurs 
according to the standard Metropolis rule. The results for all values of the magnetic 
field parameters r and a T show that finite-size effects are less-pronounced for L > 10. We 
have used 2 x 10 6 MC steps to equilibrate the system and 10 6 MC steps for averaging. 
Notice that these averages are MC or time averages. Nonetheless, in systems with 
quenched disorder we need to average over random- field realizations [23], i.e., a certain 
quantity A has the mean value given by, 

M = -£^' ( 6 ) 

1=1 

where A\ stands for MC average of the z-th sample and n r is the number of random-field 
realizations. Thus, in addition to the MC steps used for time averaging, we have used 
number of samples than went up to 2000 samples of the distributions Pi jS (H), although 
we afterwards statistically verified by means of the resampling procedure [21] that this 
number can be substantially diminished without jeopardizing the results regarding the 
analysis of the critical nature and frontiers. However, the reader must be aware that for 
a survey on the universality class of the problem, which is defined by the set of critical 
exponents, a whopping number of samples is mandatory. 

To obtain the quenched random fields we have used standard procedures for 
generating t-Student and r-distribution variables [25] . 

III. Numerical Analysis 

As previously discussed [H], the PDFs of Eqs. fl2J and (J3J may be platykurtic, for r < 1, 
or leptokurtic, for r > 1. To study the critical properties of the model we performed 
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Figure 2: (Color online) Magnetization per spin as a function of temperature for r = 2, 
<7 T = 1.2, L = 16 and some blocks with 200 samples each one. In the inset we show 
the corresponding average magnetization per spin, considering all the 2000 samples. 



MC simulations for r=— oo, 0.2, 3/2, 5/3 and 2.0, in the range 0.05 < T < 5.0 and 
0.1 < <y T < 2.4. Notice that for r > 5/3, we must use the PDF width a T instead of the 
standard deviation a. Next, we will present our results for the order parameter (m), 
the magnetization per spin, and the magnetic susceptibility \- The latter is obtained 
from the simulations by making use of the fluctuation- dissipation relation, 



where k is the Boltzmann constant (we set k — 1) and ( ) stands for averages given 
by Eq. (JOJ) - Before moving forwards we discuss the choice of our sampling. To that 
aim, we considered a quite broad distribution, r = 2 which is equivalent to the Lorentz 
distribution, and performed a total of 2000 samples in a system with L = 16 on which 
we have applied a jackknife-like procedure, which belongs to the established class of 
resampling methods in Statistics [26]. Explicitly, we regrouped our magnetization data 
into smaller groups. The smallest group size we considered was 200 elements. Already 
at this group size we could achieve statistical significance in Student t-tests for a p- value 
equal to 0.05 when we appraised the likelihood of having groups of samples of equal 
and unequal size (with non-coincident elements obtained from the same r, a T and T 
parameters) with the same statistical moments. Bear in mind that if we consider 200 
samples of a L = 8 system, we are actually taking into reckoning over than 10 5 generated 
quenched random fields, which are sufficient to have a worthy description of the fat tail 
behaviour, namely the effect of very unlikely random fields. Such a validity can still be 
intuitively grasped when we plot the magnetization curve considering different group 
sizes in Fig. [2j These data are results for simulations for r = 2, L = 16 and different 
sub-samples of 200 samples each. It is evident the concurrence between the sub-sets. In 
the inset of the same figure we plot the mean magnetization averaged over all the 2000 
samples which presents the same behaviour as the smaller sets. Once more, we would 
like to emphasize that we are not setting our sights on the critical universality of the 
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Figure 3: (Color online) Magnetization per spin as a function of temperature for a 
platykurtic case corresponding to the uniform distribution with L = 16 and typical 
values of er T , showing that ferromagnetic order occurs at finite temperatures. For 
values of a T at least up to 2.4, we observe a continuous transition between the ordered 
and the disordered phases (left side). Observe the small variation of the transition 
temperatures for increasing values of a T as it happens in the mean-field limit j!4) . It 
is also shown the Binder cumulant for a T = 2.4 and different lattice sizes (right side). 
Notice the crossing of the curves at a critical temperature near T = 3.5. 



system, but in the definition of critical regions instead, as done in previous cases (see 
e.g. Ref. [27]). 

As made in the mean-field case [H], we will study the platykurtic (r < 1) and the 
leptokurtic (r > 1) cases separately. 

III.l. Platykurtic case: r < 1 

We start presenting the results for the lower limiting platykurtic case, i.e., the uniform 
distribution corresponding to r — > oo (see Fig. EJ). The numerical results suggest that 
the transition temperature is not much affected by the values of a T , in agreement with 
mean-field predictions |T3] . In addition, plotting the data for the Binder cumulant (left 
panel of Fig. [3])), obtained by 

U l = 1 - 4% > (8) 

we can observe the signature of a phase transition, i.e., the crossing of the curves for 
different lattice sizes at the critical temperature. For 3D, comparing the latter results 
with mean-field estimates, we are still able to perceive a phase transition for widths of 
the distribution greater than 1, which is the maximum value of a T wherein symmetry- 
breaking exists in the uniformly distributed external random-field case. 

In Fig. HJ we present our results for the magnetization per spin versus the 
temperature for r = 0.2, L — 16 and typical values of a T . One can observe that 
ferromagnetic order occurs at low temperatures. We have continuous phase transitions 
between the ordered ferromagnetic phase (m — > 1) and the disordered paramagnetic one 
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Figure 4: (Color online) Magnetization as a function of temperature for another 
platykurtic case (r = 0.2) and typical values of a T (left side). It is also shown the 
Binder cumulant for tr T = 1.2 and several lattice sizes (left panel), where we can see 
that the critical temperature is located somewhere near T = 4.1. 

(m y 0) for all studied values of a T . In the right panel of the same figure, we depict 
the scale-behaviour of the Binder cumulant for a particular value of width, a T = 1.2. 
Collating both platykurtic cases we observe that the critical temperature wanes as the 
distribution approaches the Gaussian which is understandable since by increasing r we 
soar the disorder. As mentioned for the uniform case, the cubic lattice presents a larger 
interval of values of a T in which we can verify a phase transition. 

In conclusion, if the random fields are generated by platykurtic distributions 
(r < 1), the system undergoes a phase transition for all studied values of cr T , 
with a ferromagnetic phase at low temperatures and a paramagnetic phase at high 
temperatures. For large values of a T , the system is in the paramagnetic phase for all 
temperatures. In addition, the critical temperatures are not much affected by increasing 
values of r. 

III. 2. Leptokurtic case: r > 1 

The picture remains very much the same when we study systems presenting fat tailed 
behaviour in the probability density function of the quenched random fields. Namely, we 
have found evidence for a continuous phase transition for all the values of the parameter 
t we simulated. However, the critical temperatures are more sensitive to variations 
in <j T than the platykurtic region, in agreement with mean-field predictions |14j . 
Concomitantly, we have understood a decrease of the maximum value of <7 r , <r™ ax , for 
which it is possible to identify the existence of a ferro-paramagnetic transition at finite 
temperature. Once again, it is straightforward to verify that the greater r, the smaller 
cr™ ax , since an augment of the value of r encloses an increase of the disorder. For 
instance, comparing the results for r = 5/3 and r = 2, we notice that fixing a T = 2.4, 
we distinguish a finite temperature phase transition for the former r while we cannot 
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Figure 5: (Color online) Magnetization as a function of temperature for the leptokurtic 
case (r = 5/3), L = 16 and typical values of a T . The results suggest that for a T > 2.4 
the system is in the paramagnetic phase. 



descry the existence of critical behaviour at finite temperature for the latter (see Figs. E] 
and[6j). Comparing with high-dimensional result (mean field), we have the limit width 
at T = for 0>=5/3 = l/y/2 and a T= 2 = 2/n [14]. We must note that the r = 2 
corresponds to the sub-domain of values of the tail index we obtain a divergent free 
energy at T = in the mean-field limit was calculated. 

One of the points that called our attention is the fact that for large values of the 
disorder, particularly when there is an obvious tail effect, the magnetization we found 
at low temperatures is significantly different from 1. In order to check this result we 
made use of the mean-field treatment. Explicitly, we compared the numerical results 
from our Monte Carlo simulations with the solutions to the mean-field equation [T4] . 



The results of two cases are depicted in Fig. [JJ For a very small value of a T (0.1), 
the low-temperature magnetization per spin is equal to 1, as in the three-dimensional 
case, but it is visible that for a T = 0.4 we can have m < 0.9 at low temperatures. The 
integration was carried out using the global adaptive strategy algorithm [25], which for 
the two cases exhibited has an maximum error of 0.4% in the numerical integration. 
Moreover, for large values of the temperature, which imply a paramagnetic phase, we 
obtained values of |m| = 0.003, which swayed us to consider that this effect as a signature 
of the system, particularly of the tail. 

With these and other results we obtained, it is possible to sketched the diagram of 
the temperature T versus r for a fixed value o T = 1.2 (see Fig. [HJ. The squares are MC 
estimates of the critical temperatures, and we have considered the critical temperature 
of references [291 EZ] for the Gaussian case (r = 1). The computational results are in 
good agreement with the heuristic expression, 
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Figure 6: (Color online) Magnetization as a function of temperature for another 
leptokurtic case case (r = 2) a T = 1.2 and several values of the system size (left 
panel); Magnetic susceptibility as a function of temperature for another leptokurtic 
case case (r = 2) a T = 0.4 and several values of the system size (right panel). 




Figure 7: (Color online) Mean-field magnetization as a function of temperature for 
r = 2 and o-i = 0.1 and 0.4. We can observe in this high-dimensional limit the same 
behaviour at low temperatures presented by the three-dimensional case. The points 
were obtained numerically and the line is presented as guide to the eye. 

when a = 5.08 ± 0.16, b = 0.44 ± 0.04 and r c = 2.004 ± 0.002 (R 2 = 0.998 and 
X 2 /n = 0.007). This curve is shown in Fig. [8] as a dashed line. Observe that for r > 1 
the critical temperatures decrease rapidly with increasing values of r, but they do not 
decrease significantly for r < 1, as above-discussed. The validity of the adjustment 
can be tested by comparing the asymptotic limit of our ansatz (flOi) . which is equal 
to a, with the critical temperature of the uniform distribution case, T c = 4.41. This 
approximately corresponds to a 10% difference (considering error) which is acceptable 
taking into account the reduced number of parameters we assumed so that the Akaike 
information criterion is implicitly minimized [30j . 
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Figure 8: (Color online) Sketch of the phase diagram of the model in the 
plane Temperature versus r, for a T — 1.2, separating the Ferromagnetic (F) and 
Paramagnetic (P) phases. The points are Monte Carlo estimates of the critical 
temperatures and and the line is the fit using Eq. (|10[) , 
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IV. Remarks and Outlook 

In this work, we have presented a preliminary survey on both the existence and the 
nature of critical behaviour of the random-field Ising model on cubic lattices with 
nearest-neighbour interactions by means of Monte Carlo simulations based on the 
standard Metropolis algorithm. In order to generate the quenched random fields, 
we have considered a family of continuous probability density functions, defined by a 
parameter r comprising the r-distribution, for r < 1, and the Student-t, for r > 1. This 
set of distributions allowed us to study a larger variety of cases that go from the uniform 
distribution to fat-tailed distributions with the same asymptotic regime as the Levy 
distribution that were already used in the study of vitreous systems [15] , and for which 
it was recently found that the local fields follow non-Gaussian distributions [T6| IP?] . 
These models are considered relevant to the description of spin-glasses with Ruderman- 
Kittel-Kasuya-Yosida interaction and spin glasses exhibiting a wide hierarchy of coupling 
intensities as well. The current results show that in the cubic lattice, this RFIM model 
presents many of the traits found for the mean-field case [T3j. Namely, based on the 
analysis of different quantities such as the order parameter, response function and Binder 
cumulant, we have found evidence for the existence of continuous phase transitions for 
all values of r with the critical temperature, T c , depending on the characteristic width 
of the distribution, a T , and the tail index, r, itself. As expected, the spell of a T in which 
we can detect a ferro-paramagnetic phase transition is greater than the region found in 
the mean-field approach. Although we might be tempted to think that the parameter 
r represents just a supplementary way of tuning the non-thermal disorder introduced 
in the system, we must take into account results obtained at T = where the thermal 
disorder is quelled. Our previous calculations have shown quite different behaviour of 
thermodynamical quantities depending on the parameter r, particularly for distributions 
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fatter than the Lorentzian. We expect that the analysis of the ground state together 
with upgraded simulations we are also carrying out manage to shed further light on the 
actual impact of the tails. 
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